LIMITS OF BIFRACTIONAL BROWNIAN NOISES 



MAKOTO MAEJIMA AND CIPRIAN A. TUDOR 

Abstract. Let B H - K = (^B^' K ,t > 0^ be a bifractional Brownian motion 
with two parameters H € (0, 1) and K £ (0, 1]. The main result of this paper 
is that the increment process generated by the bifractional Brownian motion 
( B K+t - B h' K ' l ^ °) converges when h -> oo to (2( 1 - K )/ 2 Bf K , t > 0), 
where (B^ K t t > 0) is the fractional Brownian motion with Hurst index HK. 
We also study the behavior of the noise associated to the bifractional Brow- 
nian motion and limit theorems to B H ' K . 



1. Introduction 

Introduced in [4], the bifractional Brownian motion, a generalization of the 
fractional Brownian motion, has been studied in many aspects (see [1], [3], [6], 
[7], [8], [9] and [10]). This stochastic process is defined as follows. Let H G (0, 1) 

and K G (0, 1]. Then B H - K = (B?' K ,t> o) is a centered Gaussian process with 
covariance 



E 



B" ,K B H ' K 



= 2- K ((t 2H + s 2H ) K -\t- s\ 2HK ) 



When K = 1, it is the fractional Brownian motion B H = [B^ ,t > 0) with the 
Hurst index H G (0, 1). In general, the process B H ' K has the following basic prop- 
erties: It is a selfsimilar stochastic process of order HK £ (0, 1), the increments 
are not stationary and it is a quasi- helix in the sense of [5] since for every s,t > 0, 
we have 



2- K \t-s\ 2HK <E 



K _ gH,K 



<2 1 - K \t-s\ 2HK . 



The trajectories of the process B H ' K are 5-Holder continuous for any 6 < HK and 
they are nowhere differentiable. 

A better understanding of this process has been presented in the paper [7], 
where the authors showed a decomposition of B H ' K with H, K G (0, 1) as follows. 
Let (Wg, 6 > 0) be a standard Brownian motion independent of B H,K . For any 
K G (0, 1), they defined a centered Gaussian process X K = (X^ , t > 0) by 

X? = \ (1 - e-*<)0- (1+K)/ W 9 . (1.1) 
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Its covariance is 

E [X t K Xf] = T(l - K)K- X (t K + s K -(t + s) K ) . (1.2) 
Then they showed, by setting 

X^ K :=X« H , (1.3) 

that 

(CiX?'* + Bf t > 0) 4 (C 2 Bf K , t > 0) , (1.4) 

where Ci = (2- x if(r(l - K))- 1 ) 1 / 2 , C 2 = 2^- K ^ 2 and = means equality of all 
finite dimensional distributions. 

The main purpose of this paper is to study the increment process 

(BZ$-B»' K ,t>0) 

(where h > 0) of B H ' K and the noise generated by B H ' K and to see how close this 
process is to a process with stationary increments. In principle, since the bifrac- 
tional Brownian motion is not a process with stationary increments, its increment 
process depends on h. But in this paper we show, by using the decomposition 
(1.4), that for h — > oo the increment process (^B^^ — B^' K ,t > 0^ converges, 
modulo a constant, to the fractional Brownian motion with Hurst index HK in 
the sense of finite dimensional distributions, so the dependence of the increment 
process depending on h decreases for very large h. Somehow, one can interpret 
that, for very big starting point, the bifractional Brownian motion has stationary 
increments. Then we will try to understand this property from the perspective of 
the "noise" generated by B H ' K i.e. the Gaussian sequence B^^ — B^' K , where 
n > are integers. The behavior of the sequence 



Y a (n) = E [(Sf'f - flf (B»f n+1 - B»f n ) 



a e 



(which, if K = 1, is constant with respect to a and of order n 2H 2 ) is studied 
with respect to a and with respect to n in order to understand the contributions 
of B HK and X H > K . 

We organize our paper as follows. In Section 2 we prove our principal result 
which says that the increment process of B H ' K converges to the fractional Brow- 
nian motion B HK . Sections 3-5 contain some consequences and different views 
of this main result. We analyze the noise generated by the bifractional Brownian 
motion and we study its asymptotic behavior and we interpret the process X H ' K 
as the difference between "the even part" and "the odd part" of the fractional 
Brownian motion. Finally, in Section 6 we prove limit theorems to the bifractional 
Brownian from a correlated non-stationary Gaussian sequence. 



2. The limiting process of the bifractional Brownian motion 

In this section, we prove the following main result; it says that the increment 
process of the bifractional Brownian motion converges to the fractional Brownian 
motion with Hurst index HK. 
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Theorem 2.1. Let K e (0, 1). Then 

(<+f - B»> K ,t > 0) 4 ^-^Bf K ,t > 

where =4> means convergence of all finite dimensional distributions. 

To prove Theorem 2.1, we use the decomposition (1.4). It is enough to show 
that the increment process associated to X HK converges to zero; we prove it in 
the next result, and actually we measure how fast it tends to zero with respect to 
L 2 norm. It will be useful to compare this rate of convergence with results in the 
following sections. 

Proposition 2.2. Let X H,K be the process given by (1.3). Then as h — > oo 

r(l - K)K- 1 2 K H 2 K(l - K)t 2 h 2(HK - 1 \l + o(l)). 



E 



f Y H,K vH,K\ 



As a consequence, 



X 



h+t 



X*' K ,t> 



(X(t) = 0,t>0) ash^oo. 



Proof. Note from (1.2) and (1.3) that 



E 



X HK X H ' K 



r(i - k)k- x (t 2HK + s 



2HK 



(t 2H + s 2H 



and in particular, for every t > 

21 



We have 



E 



X 



H.K 



X 



E 



H.K 



E 



r(l- K)K- 1 {2-2 K )t 2HK 



2E 



V H,K Y H,K 



E 



(x^-xry 



(2 - 2 K )h 2HK ) 



L /i+t 

Then 

/ :=K(T(1 -K))-^ 

= ({2-2 K )(h + t) 2HK 

- 2 ({h + t) 2HK + h 2HK - {{h + t) 2H + h 2H ) K ) + 

= -2 K ((h + t) 2HK + h 2HK ) +2((h + t) 2H + h 2H ) K 

= -2 K h 2HK (1 + (th- 1 ) 2HK + 1) + 2h 2HK ((1 + th- x ) 2H + iy 
Therefore for very large h > we obtain by using Taylor's expansion 
/ = -2 K h 2HK (2 + 2HKth- 1 + H(2H - l)t 2 /i~ 2 (l + o(l))) 

+ 2h 2HK (2 + 211th- 1 + H(2H - l)t 2 h- 2 (l + o(l)))' 

Now we use Taylor expansion for the function (2 + Z) K for Z close to zero. In 
our case Z = 2Hth~ 1 + H(2H - \)t 2 hr 2 + r(h) with r(h)h 2 -» as h -> oo. We 
obtain 

/ = -2 K h 2HK (2 + 2HKth- 1 + H(2H - 1)< 2 ^" 2 (1 + o(l))) 



, K 
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+ 2h 2HK (2 K + K2 K - 1 (2Hth- 1 + H(2H - l)t 2 h- 2 + r(h)) 

+ 2- 1 K(K - l)2 K - 2 (2Hth- 1 + H(2H - l)* 2 ^ 2 + r{h)f + o{hr 2 )^ 
= h 2HK 2 K HK(-2HK + 1 + 2H - 1 + UK - H)t 2 h- 2 {1 + o(l)) 
= h 2HK 2 K H 2 K(l - K)t 2 hr 2 {\ + o(l)). 
Consequently, we have 



E 



r(l - K)K- x 2 K n 2 K{\ - K)t 2 h 2i - HK - 1 \l + o(l)), 



which tends to as h — > cxo, since Ifif — 1 < 0. □ 



3. Bifractional Brownian noise 



By considering the bifractional Brownian noise, which are increments of bifrac- 
tional Brownian motion, we can understand Theorem 2.1 in a different way. Define 
for every integer n > 0, the bifractional Brownian noise 



Y„ 



R H,K pH,K 
B n+l ~ B n ■ 



Remark 3.1. Recall that in the fractional Brownian motion case (K = 1) we have 
for every sgM and for every n > 0, E [Y a Y a+n ] = E [Y Y n ] . 



Let us denote 



R(0, n) = E[Y Y n ] = E [b™ (b%* B^ k ) 



and 



R(a, a + n)=E [Y a Y a+n ] = E [(b»£ B^ k ) (b»£ +1 B?£)] . (3.1) 

Let us compute the term R(a, a + n) and understand how different it is from 
R(0, n). We have for every n > 1, 



\2H 



K 



— n 



2HK 



R(a, a + n) = 2~ K (((a + l) 2H + {a + n + l) 2 

- ((a + l) 2H + (a + n) 2H )* - (n - l) 2HK 

- (a 2H + (a + n+ l) 2H ) K - (n + l) 2HK 
+ (a 2H + {a + n) 2H ) 1 

=:2- K (f a (n)+g(n)), 



ls n 2HK ) 



(3.2) 



where 



j a (n)= {{a+l) 2H + (a + n+lf 



((a + l) 2H + (a + n) 2H y 



(a 2H + (a + n+ l) 2H ) K + (a 2H + (a + n) 2H ) 



K 
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and for every n > 1 , 



g{n) = {n+l) ztlK + (n - I) 



\2HK 



- 2n 



2HK 



Remark 3.2. (i) The function g is, modulo a constant, the covariancc function of 
the fractional Brownian noise with Hurst index HK. Indeed, for n > 1, 



g(n) = 2E [B? K (B™-B* 



}HK\ 



(3.3) 



(ii) g vanishes if 2HK = 1. 

(iii) The function f a is a "new function" specific to the bifractional Brownian case. 
(Note that f a vanishes in the case K = 1.) It corresponds to the noise generated 
by X H ' K . Indeed, it follows easily from (1.4) that 



/„(n) = -2 K C 2 E 

=: -2 K C\R X 
for every a and n E N. 



/ V H,K vH,K\ ( vH.K vH,K\ 

(a, a + n) 



(3.4) 



We need to analyze the function f a to understand "how far" the bifractional 
Brownian noise is from the fractional Brownian noise. In other words, how far is 
the bifractional Brownian motion from a process with stationary increments? 

Theorem 3.3. For each n it holds that as a — > oo 

f a (n) = 2H 2 K(K - l)a^ HK - 1 Hl + o(l)). 
Therefore lim f a (ri) — for each n. 

The bifractional Brownian noise is not stationary. However, the meaning of the 
theorem above is that it converges to a stationary sequence. 



Proof. We have, for a — > oo, 



f a (n) = a 2HK {(1 + a^f H + (1 + („ + lja" 1 ) 2 "}* 



{(1 + a- 1 ) 2 * + (1 + m- 1 )^}^ -{! + (! + („+ l)a~ 1 ) 2H } K 



= a 



+ {1 + (1 + 

2HK 



[{1 + 2Ha~ 1 + H(2H - l)a- 2 {l + o(l)) 

+ 1 + 2if(n + lja" 1 + H{2H - l)(n+ l) 2 a~ 2 (l + o(l))} K 

- {1 + 2 J ffa" 1 + if(2if - l)a~ 2 (l + o(l)) 

+ 1 + 2Hna- 1 + H{2H - l)n 2 a~ 2 (l + o(l))} K 
-{1 + 1 + 2H{n + l)a _1 + H(2H - l)(n + l) 2 a~ 2 (l + o(l))} K 

- {1 + 1 + 2ff (n + l)a _1 + H(2H - l)(n + l) 2 a~ 2 (l + o(l))} K 



2a 



2HK 



{1 + H(n + 2)a- 



MAKOTO MAE.TIMA AND CIPRIAN A. TUDOR 



+ 2~ 1 H(2H - 1)(1 + (n + l) 2 )a" 2 (l + o(l))} K 

- {l + if(n + l)a _1 + 2- 1 H(2H - 1)(1 + n 2 )a~ 2 (l + o(l))} 

- {1 + ff(n + lja" 1 +2~ 1 H(2H - l)(n+ l) 2 a~ 2 (l + o(l))} 



A" 



A" 



1 + HnaT 1 + 2~ 1 H(2H - l)n 2 a" 2 (l + o(l))} 



A 



2a 



2HK 



{l + K(H(n + 2)a- 1 

+ 2- 1 H(2H - 1)(1 + (n + l) 2 )a- 2 (l + o(l))) 

+ 2"^ - l)(H(n + 2)a- 1 (l + o(l))) 2 (l + o(l))} 

- {1 + /T(i2"(n + l)a _1 + 2~ 1 H(2H - 1)(1 + n 2 )a~ 2 (l + o(l))) 

+ 2- 1 tf(A' - l)(#(n + ljo-^l + o(l))) 2 (l + o(l))} 

- {l + if(ff(n+ l)^ 1 

+ 2~ 1 H{2H - 1)(1 + (n + l) 2 )a" 2 (l + o(l))) 
+ 2"^ (# - l)(#(n + ljo-^l + o(l))) 2 (l + o(l))} 
+ {1 + KiHna- 1 + 2~ 1 H(2H - l)n 2 a~ 2 (l + o(l))) 

+ 2~ 1 K(K - l^Hna-^l + o(l))) 2 (l + o(l))}" 



= 2a 2HK (if#(n + 2) - Kff(n + 1) - Ki?(n + 1) + KHn)^ 1 

+ 2~ 1 KH(2H - 1)(1 + (n + l) 2 ) + 2~ 1 K(K - l)H 2 {\ + (n+ l) 2 ) 

- 2- 1 KH(2H - l)(n 2 + 1) - 2- 1 K(K - l)H 2 (n + l) 2 

- 2- 1 KH(2H - l)(n + l) 2 - 2- 1 K(K - l)H 2 (n + l) 2 

+ 2~ 1 KH(2H - l)n 2 + 2~ 1 K{K - l)H 2 n 2 )} a - 2 (l + o(l)) 

= 2H 2 K(K - l)a 2{HK - 1 \\ + o(l)). 
Since iJX — 1 < 0, the last term tends to when a goes to the infinity. □ 



Remark 3.4. The fact that the term / a (n) converges to zero as a -> oo could be 
seen by Proposition 2.2 since, using Holder inequalities, 



R xH ' K (a, a + n) < ( E 



K vH.K 

+ 1 ~ A a ' 



1/2 



/ X H,K 



x. 



H.K 
a+n 



1/2 



and both factors on the right hand side above are of order a HK 1 . The result 
actually confirms that for large a, X^_^ +1 — X^^ is very close to X^^ — X^ ,K ■ 



4. The behavior of increments of the bifractional Brownian motion 

In this section we continue the study of the bifractional Brownian noise (3.1). 
We are now interested in the behavior with respect to n (as n — > oo). We know 



BIFRACTIONAL NOISES 



7 



that as n — > oo the fractional Brownian noise with Hurst index HK behaves 
as HK(2HK — l)n 2 ^ HK ~ 1 h Given the decomposition (1.4) it is natural to ask 
what the contribution of the bifractional Brownian noise to this is and what the 
contribution of the process X H K is. We have the following. 

Theorem 4.1. For integers a, n > 0, let R(a, a + n) be given by (3.1). Then for 
large n, 

R(a, a + n) = 2~ K (2HK(2HK - l)n 2(HK -^ 

+HK(K - 1) ((a + 1) 2H - a 2H ) n 2(^-i)+(i-2H) + . . . y 

Proof. Recall first that by (3.2) and (3.3), 

R(a,a + n) = 2- K (f a (n)+g(n)) 

and the term g(n) behaves as 2H K (2H K — l)n 2 ( HK ~ 1 ) for large n. Let us study 
the behavior of the term f a (n) for large n. We have 

f a {n) = ((a + l) 2H + (a + n+ l) 2H f ((a + l) 2H + (a + n) 2H f 



(a 2H + (a + n+ l) 2H ) K + (a 2H + (a + n) 2H ) 

2H\ 1 



n 



2HK 



K 



(((a+l)n- 1 ) 2H + ((o + ljn^ + l) 

- ({{a+l)n- l ) 2H + {an- 1 + l) 2H ) 

- ((an- 1 ) 2 " + ((a + l)n- 1 + l) 2H ) 
+ ((an- 1 ) 2H + (an- 1 + l) 2H ) K 



K 



K 



2HK 



2H-2H 



(a + l) 2H n 



+ 1 



+ iyr^H(2H -1)...(2H- j)(a + 1)' +1 ? 

3=0 



-3-1 



K 



- [{a + l) 2H n - 2H + 1 

oo K 

+ 2(0' + iVr^HiZH -1)...(2H- jy^n-i- 1 ) 

4— n / 
+ 1 

OO \ 

]T(0 + iy-r^H(2H -1)...(2H- j)(a + l^n"'" 1 J 

4— n / 
\ i (r",r-" + 1 



3=0 
2H-2H 



— a n 



K 
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+ £(c? + i v-y l2H ( 2H - 1) ■ ■ ■ (2^ - 

J=0 



i+i n -j-i 



By the asymptotic behavior of the function (1 + y) K when i/^Owe obtain 



Ja(n) = n 



2HK 



oo 



l + ^I + l)!)- 1 ^-!)-^^ 

t=o 

> ( (a+l) 2H n- 2H + 1 

OO \ 

+ 1)\)- 1 2H(2H -1)...(2H- j)(a + l^n"'" 1 J 



3=0 



£=0 



x ( (a+l) 2H n- 2H + l 



+ l)!)"^ (2JT - l) • • • (2H - 3 )a^ 1 n-^ 1 

3=0 



t+1 



-l-^l + l)l)-^K(K-l)...(K-£) 



x | a 2H n- 2ff + 1 



OO \ 

+ + 1)\)~ 1 2H(2H -1)...(2H- 3 )a^n-i- x 

3=0 ' 



t+1 



1=0 



x I a 2H n- 2H + 1 



OO s 

+ X)(C? + iy.)- 1 2H(2H -1)...(2H- j)a? + *n-t- x ) 



t+i 



= n 



2HK 



1 + J2((l + m- 1 K(K -!)■■■ (K 



x (a+l) 2H n- 2H +1 
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oo v e+i 

+ + 1)\)- 1 2H(2H -1)...{2H- j)(a + l^n"'" 1 J 

j=0 ' 

-l-Y j {{l + l)y)-'K{K-l)...{K-t) 
x [{a+lf H n- 2H + 1 

1X3 s £+1 

+ ]T((j + l)!)- 1 2^(2if - 1) • • • (2H - jy+'n-i- 1 J 

oo 

e+i 

+ + l)!)"^ (2*T - 1) • • • (2ff - j)^ 1 ^" 1 

3=0 



l + 1 £ i ((l + l)\)- 1 K(K-l)...(K 



<=1 



x I a 2ff n- 2ff + 1 



oo s e+i 

+ l)!)- 1 2ff(2if - 1) • ■ • (2ff - jJ^'+^-J- 1 J 



= 2" 1 J ftT(X - l)n 2HK [ ((a + l) 2ff n- 2H + 2ff (a + I)™" 1 

+H(2H -l)(a + l) 2 n- 2 ) 2 

- ((a + l) 2H n- 2H + 2Han- 1 + H(2H - l)a 2 n - 2 ) 2 

- (a 2H n- 2H + H(a + l)^ 1 + 2H(2H - l)(a + I) 2 "" 2 ) 
+ (a 2H n- 2H + Han- 1 + 2H(2H - l)a 2 n - 2 ) 
+ ••• 

= HK(K - 1) ((a + l) 2ff - a 2if ) n ^HK-i)+{^H) + . . . _ 
This completes the proof. □ 



Let us discuss some consequences of the theorem above. 

Remark 4.2. What is the main term in R(a,a + n)? Note that H > ^ if and 
only if 2(HK - 1) > 2(#if - 1) + (1 - 2if). Consequently the dominant term for 
R(a, a + n) is of order n 2HK - 2 if H > ± and of order n ^HK-i-2H \f H < \. 
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Another interesting observation is that, although the main term of the covari- 
nace function R(a, a + n) changes depending on whether H is bigger or less than 
one half, the long-range dependence of the process B H ' K depends on the value of 
the product HK. 

Corollary 4.3. For integers a > 1 and n > 0, let R(a, a + n) be given by (3.1). 
Then for every a E N, we have 

R( a , a + n) = oo if 2HK > 1 

n>0 

and 

R(a, a + n) < oo if2HK<\. 

n>0 

Proof. Suppose first that 2HK > 1. Then it forces H to be bigger than | and the 
dominant term of R(a, a + n) is n 2HK ~ 2 when n is large. So the series diverges. 

Suppose that 2HK < 1. If H > |, the main term of R(a, a + n) is n 2HK - 2 and 
the series converges. If H < ^, then the main term is n 2 #^-2ff-i an( j jjjg se ries 
converges again. 

If 2HK = 1 (and thus H > |) then R(a, a + n) behaves as n — > oo as n^ 1 ^ 2H 
and the series is convergent. □ 

Corollary 4.4. Let R xH K (a, a + n) be the noise defuned in (3.4). Then 

R x ' (a, a + n) 

= T(l - K)K- X (-AHK(K - 1) {{a + 1) 2H - a 2H ) rj 2(^-i)+(i-2H) + . . . y 

Proof. It follows from Theorem 4.1 and the fact that the covariance function of the 
fractional Brownian motion with Hurst parameter HK behaves as HK(2HK — 



5. More on the process X H K 

We will give few additional properties of the process X K defined in (1.1). Recall 
(1.2) that for every s,t > 

R xK (a, t) := E[Xfxf\ = T(l - K)K- x {t K + s K - (t + s) K ). 

Denote by B K I 2 = (B^ 2 , t G M) a fractional Brownian motion with Hurst index 
K/2 defined for alH e R and let 



e,K,2 =2 - 1 f B K,2 + B K/2y 



The processes B°> K / 2 and B e,K / 2 are respectively the odd part and the even part 
of the fractional Brownian motion B K I 2 . Denote by R° ,K / 2 the covariance of 
the process B°' K ^ 2 , by R e,K / 2 the covariance of the B e ' K / 2 and by R B ' the 
covariance of the fractional Brownian motion B K I 2 . We have the following facts: 

R* K {f, s) = C 3 R BK/2 (t, -s) = C 3 R BK/2 (s, -t) 
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where C 3 = 2I\1 - K)^ 1 , and 

R°> K '\t,s) = \ (R BK/2 (t,s) - R BK/2 (t,-s)) 

and 

R e > K '\t,s) = \ (R BK/2 (t,s) + R BK/2 (t,- S )) . 

As a consequence 

R*> K ' 2 {t, s) - R°> K ' 2 {t, s) = R BK/2 (t, -s) = C^R xK (t 7 s). 

From the above computations, we obtain 

Proposition 5.1. We have the following equality 

C- 1/2 X K + B e ' K / 2 = B°' K ^ 2 

if X K and B e > K / 2 are independent. 

Let us go back to the bifractional Brownian noise R(a,a + n) given in (3.1). By 
the decomposition (1.4), we have 

C 1 X H ' K + B H ' K ±C 2 B HK , 



where C\ and Ci are as before, and thus we get 

R(a, a + n) = C\R bHK (a, a + n) - ClR x "' K (a, a + n) 

= C 2 R B " K (0, n) - C 3 (R e ' K / 2 ' H (a, a + n)- R°> K ' 2 > H {a, a + n)) 

where R e ' K / 2 ' H (a, a + n) is the noise of the process B e t ^ 2 \t > 0. 

Remark 5.2. The fact that the covariance function R x (a, a + n) of the process 
X K I 2 converges to zero as a — > oo can be interpreted as "the covariance of the odd 
part" C 3 R B °' K ' 2 (a, a+n) and "the covariance of the even part" C3R B "' K/2 (a, a+n) 
have the same limit 2~ 1 C 2 R BK/2 (0,n) when a tends to infinity. 



6. Limit theorems to the bifractional Brownian motion 

In this section, we prove two limit theorems to the bifractional Brownian motion. 
Throughout this section, we use the following notation. Let < H < 1,0 < K < 1 
such that 2HK > 1 and let = 1,2,...) be a sequence of standard normal 

random variables. Define a function g(x,y),x > 0,y > by 

g(x,y) = 2 2 - K H 2 K(K l)(x 2H + y 2H ) K - 2 {xy) 2H -' 

+ 2 1 - K HK(2HK -l)\x- y\ 2HK ~ 2 
=■ 9i(x, y) +92{x, y), (6.1) 



for (x, y) with x ^ y and x ^ and y ^ 0. 
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Proposition 6.1. Under the notation above, assume that E[^i^j] = g(i,j). Then 

[at] \ 

&t>0U(Bf',t>0). 



-HK ' 



To prove this, we need a lemma. 
Lemma 6.2. 



f [ S g{u, v)dudv = 2- K [(t 2H + s 2H ) K -\t- s\ 2HK ] . 
Jo Jo 

Proof. It follows easily from the fact that §^{x,y) — g{x,y) for every x,y > 
and by using that 2HK > 1. 



□ 



Proof. (Proof of Proposition 6.1.) It is enough to show that 

M 

n -HK c 1 I „-HK 



[nt] 
i=l 



In ~E 

-> E[B?' K Bf> K ] = 2~ K ((t 2H + s 2H ) K -\t- s\ 2HK ) . 



We have 



>//./\ oi?,-ffj _ ((t 2H + s 2H ) K - u /s 

[ni] [ns] [ni] [ns] 



/» = n- 2HK E E mud - n- 2HK j: E 9(i, ;)■ 
i=i j=i i=i j=i 



Observe that 



z 7 

fll - - 

n n 



-,2-i? rr2 



2ff 



2ff" 



Jf-2 



+ 2 l ~ K HK(2HK - 1) 



1 _ 1 
n n 



2HK-2 



2H-1 



= 2 2 - K H 2 K(K - l) n -2H(*-2)-2(2/f-l) ( .2H +j 2H ) ^-2 ( . j) 2H-l 

+ 2 1 - K HK(2HK - l)n- 2HK+2 \i - j\ 2HK - 2 
= n 2 ^- HK ^ (2 2 - K H 2 K(K - l)(i 2H + j 2H ) k - 2 (ij) 2H - 1 
+ 2 l - K HK(2HK - l)\i - j\ 2HK - 2 ^j 



= n 2 ^- 



HK) g{i,j). 



(6.2) 



Thus, as n — > oo, 



n£] [ns] 

2HK\^\^ 2HK-2„ ( 1 J 



EE- 

i=i j=i 



[nt] [ns] 

" 2 EEff 

»=i j=i 



-'o 



g(u,v)dudv = 2~ K ((t 2H + s 2H ) K -\t-s 



n n 

\2HK\ 



= E[B?> K B?*] 
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□ 



Remark 6.3. This result seems easy to be generalized to more general Gaussian 
selfsimilar processes such that their covariance R satisfies € L 2 ((0, oo) 2 ). 

We next consider more general sequence of nonlinear functional of standard 
normal random variables. Let / be a real valued function such that f(x) does not 
vanish on a set of positive measure, E[f(£i)] — and -E-[(/(£i)) 2 ] < oo. Let H^(x) 
denotes the fc-th Hermite polynomial with highest coefficient 1. We expand / as 
follows (see e.g. [2]): 

oo 

= E c kHk(x), 

fc=l 

where c fe^' < 00 j c k = E[f(£i)Hk(£j)]- This expansion is possible under the 

assumption Ef(£i) = and £?[(/(£i)) 2 ] < oo. Assume that c\ ^ 0. Now consider 
a new sequence 

%=/(&),j = l,2,..., 

where = 1,2, ...) is the same sequence of standard normal random variables 

as before. 

Proposition 6.4. Under the same assumptions of Proposition 6.1, we have 

[nt] \ 

A( Cl B?' K ,t>0). 



-HK ' 



3 = 1 



Proof. Note that i]j = /(£,) = a£j + Y^k=2 c fc-#fc(£j)- We have 



it] 



-HK ' 



3=1 3 = 1 

By Proposition 6.1, it is enough to show that 



-HK 



[nt] oo 

EX>^&)- 

j=l fc=2 



[nt] oo 
i=l fc=2 



as n — > oo. 



We have 



[nt] oo 
j'=l fc=2 



[nt] [nt] oo oo 

= n _MK EEEE c '*w&)]- 

i=l j=l fc=2 £=2 

In general, if £ and 77 are jointly Gaussian random variables with E\£\ = E[rj\ = 0, 
E[£ 2 } = E[n 2 } = 1 and E[£rj\ = r, then 

E[H k (S)H t (ri)] = 5 kA r k k\, 
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where 

fi, k = e, 



5k,i 
Thus 



o, k^e. 



oo [nt] oo 

1=2 i,j=l-v£j C=2 

oo [nt] oo 

1=2 i,j = l;i^j 1=2 

Since for every i,j>l(i^= j) one has \g(i, j)\ < (E^f]) 1 ^ 2 [E^ 2 ]) 1 ^ 2 = 1, we get 

oo oo [nt] 

1=2 1=2 i,j=l;ijtj 

oo / oo \ / [nt] , . . N 2 

£=2 V=2 / \ i,j=l ; i^j V 

where we have used (6.2) 

Here asm oo, since X)fc2 c <?^ < 00 an< ^ > 1 we obtain that 

tn 1 ~ 2HK Y^T=2 ° 2 ^- converges to zero as n — > oo. On the other hand, with C an 
absolute positive constant and g\ and g 2 given by (6.1), 



n 



[nt] . . . x 2 / ["*] / . . n. 2 [nt] / . . x 2 N 

2 E E 4-A + E 



The first sum n~ 2 X)!™=i-i/j Si (^) i s a Riemann sum converging to the in- 
tegral Jq g 2 (x,y)dxdy. Note that this integral is finite because |<7i(a;,t/)| < 
C(xy) HK ~ 1 and the integral J J * |a; - y\ 2HK ~ 2 dxdy is finite when 2i/if > 1. 
Since n 2l ^ HK ~ v> — > we easily get 



[nt] / . . s 2 

n 2(tf*-l) n -2 V- ^0 

^ , \n n J 

as n — > oo. 

The second sum involving (72 appears in the classical fractional Brownian case 
because it is, modulo a constant, the second derivative of the covariance of the 
fractional Brownian motion with Hurst parameter HK. The convergence of 

[nt] / . . x 2 

n 2 ( HK-l) n -2 J- 92 



has been already proved in e.g. [2]. The proof is completed. □ 
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